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A simple  relationship  between  infrared  attenuation  in  the 
8-12  ym  region  and  liquid  water  content  of  the  atmosphere  is 
derived  by  using  a straightforward  moment  analysis.  A variety 
of  experimental  and  theoretical  data  is  used  to  demonstrate  the 
validity  of  this  relationship  under  limited  visibility  condi- 
tions. 
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I.  INTRODUCTION 


The  prediction  of  multi spectral  aerosol  attenuation  in 
limited  atmospheric  visibility  as  typified  by  water  hazes,  fogs, 
and  clouds  remains  an  outstanding  problem.  In  principle  it  is 
possible,  under  well-controlled  circumstances,  to  measure  the 
optical  composition  and  size  distribution  of  atmospheric  partic- 
ulates and  then  to  use  this  information  in  a Mie  calculation  to 
predict  the  multispectral  transmission  characteristics.  In 
practice,  however,  such  a particle  size  distribution  measure- 
ment is  far  from  routine  and  proves  to  be  both  difficult  and 
inadequate  in  application  due  to  the  spatial  and  temporal 
fluctuations  associated  with  the  real  atmosphere.  It  would 
therefore  appear  to  be  more  practical,  if  not  necessary,  to 
adopt  a more  phenomenological,  albeit  approximate,  approach  to 
aerosol  modeling  that  does  not  depend  upon  direct  measurement 
or  detailed  description  of  any  particular  distribution. 

Most  current  aerosol  models  erroneously  rely  upon  a single 
representative  particle  size  distribution.  The  current  LOWTRAN 
aerosol  model  (Ref.  1),  for  example,  uses  measured  optical  prop- 
erties (representative  of  average  continental,  rural,  urban,  or 
maritime  conditions)  with  a prototypical  distribution  to  con- 
struct by  means  of  a Mie  computation  a scaling  model  for  extrap- 
olation of  transmission  in  the  visible  wavelength  range  to 
transmission  in  the  infrared.  The  underlying  assiimption  is  that 
for  a particular  atmospheric  condition,  such  as  a continental 
water  haze,  the  shape  or  functional  form  is  a fixed  and  well- 
known  distribution.  In  many,  if  not  most,  cases  this  is  not  a 
valid  representation.  Unfortunately,  the  shape  or  functional 


form  of  the  distribution  varies  over  such  a wile  range  that  a 
prototypical  distribution,  such  a'-  a continental  haze,  is  not 
very  representative  of  aerosol  conditions  at  any  particular 
time  or  place.  For  example,  in  an  evolving  fog  formation  the 
water  droplet  distribution  tends  to  grow  in  the  sense  that  there 
are  relatively  more  large  particles  as  the  visibility  becomes 
lower.  This  is  illustrated  dramatically  in  Pig.  1,  where  we 
have  plotted  some  representative  measurements  made  in 
Grafenwohr,  Federal  Republic  of  Germany,  over  a 12-hour  period 
on  30-31  December  1975  (Ref.  2).  Figure  2 shov/s  the  same  prin- 
ciple for  particle  size  distributions  for  hazes  and  radiation 
fogs  from  the  recent  review  of  Tomasi  and  Tampieri  (Ref.  3). 

Each  of  these  distributions  leads  to  a different  spectral  de- 
pendence for  the  aerosol  extinction  coefficient,  as  shown  ex- 
plicitly in  Fig.  1. 

In  this  paper  we  propose  a simple  but  general  scaling 
model  that  takes  into  account  such  changes  for  8-12  pm  infrared 
extinction  but  does  not  depend  explicitly  upon  the  selection  of 
any  single  measured  or  assumed  size  distribution.  Although  the 
application  in  this  paper  will  be  primarily  to  water  droplet  or 
fog/cloud  optical  propvrties,  the  arguments  used  to  postulate 
she  scaling  model  are  generally  valid  also  for  other  particu- 
lates . 
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PARTICLE  RADIUS,  MICRONS 


FIGURE  i.  Growth  of  aerosol  droplets  as  fog  builds  over  time, 
Grafenwohr,  FRG,  30-31  December  1975.  The  corre- 
sponding change  in  measured  transmission  is  noted 
accordingly.  The  visibility  as  well  as  the  scaling 
between  different  spectral  bands  changes  dramatical- 
ly as  the  particle  size  dirtribution  evolves,  with  a 
corresponding  increase  in  larger  particles.  Since 
the  actual  particle  distribution  measurements  sample 
a particle  count  rather  than  volume,  we  have  plotted 
n(r)  rather  than  r3n(r)  in  Figs.  1 and  2.  A similar 
plot  of  r3n(r)  would  show  maxima  towards  larger  par- 
ticles than  n(r).  For  example,  the  case  at  10:00 
hours  ( ) would  show  a maximum  at  about  4 ym. 


LOG  OF  PARTICLE  RADIUS.  MICRONS 


FIGURE  2 


Particle  size  distributions  representative  of  hazes 
and  radiation  fogs  (Ref.  3). 


■1 

■i 


II. 


DISCUSSION 


5 


i 


To  develop  a more  flexible  model  we  will  recast  the 
aerosol  extinction  theory  in  terms  of  a moment  analysis.  If 
one  can  measure  or  predict  both  the  entire  particle  size  dis- 
tribution n(r),  where  r is  the  particle  radius,  and  the  optical 
properties  of  the  particulates  as  characterized  by  the  complex 
index  of  refraction,  then  ideally  the  problem  of  predicting 
the  multispectral  aerosol  extinction  coefficients  is  solved  by 
an  application  of  Mie  theory.  Alternatively,  one  can  describe 
the  entire  distribution  by  an  equivalent  knowledge  of  the  com- 

/If 

r r.(r)dr  for  k = 0,  1,  2 ....  One 
can  thus  recast  the  problem  by  stating  that  if  the  set  of 
moments  for  the  distribution  as  well  as  optical  proparties  are 
given,  the  extinction  coefficients  can  also  be  obtained.  For 
this  equivalence  to  have  any  oracticai  significance  we  m.ust 
show  that  the  set  of  moments,  or  preferably  a limited  selection 
of  moments,  plus  optical  properties  can  provide  useful  infor- 
mation concerning  the  aerosol  extinction  coefficient  B . We 
would  like  to  show  in  particular  that  one  specified  moment, 
namely  the  third  moment  or  the  liquid  water  content  or 
volume  density  v,  is  especially  useful  for  characterizing  in- 
frared propagation. 


For  a specified  wavelength  X that  is  large  compared  to 
the  typical  particle  radius  r,  the  Kie  extinction  cross  section 
is  proportional  to  r"*  ~ r-*)  in  an  absorptive  medium  (Ref, 


The  total  aerosol  extinction  coefficient  6 is  given  by 

aer 


the  average  of  the  cross  section  over  the  particle  size  distri- 
bution n(r).  Hence,  for  X large  compared  to  r,  which  should 
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be  appropriate  for  far  IR  radiation  and  small-particle  fogs  or 
hazes,  the  aerosol  attenuation  at  a given  wavelength  depends 
primarily  upon  the  third  moment  of  the  distribution,  or 
specifically  th3  total  volume  or  liquid  content  v of  particu- 
lates, in  a straightforward  linear  fashion:  ^3  ~ 

The  constant  of  proportionality  is  dependent  upon  the  explicit 
wavelength  due  to  its  functional  relationship  to  the  complex 
index  of  refraction.  If  we  next  examine  the  other  extreme, 
sometimes  referred  to  as  the  geometric  limit,  where  the  parti- 
cles are  typically  large  compared  to  the  wavelength  of  the 

radiation,  then  large-particle  scattering  dominates*  (Ref.  ^), 

2 

and  we  have  = 2irr  . Once  again,  a size  averaging  yields 

the  desired  result,  related  now  to  the  effective  area  A or 
second  moment  I2  for  the  given  distribution:  ~ 

Since  there  is  such  a small  change  (from  to  I2)  in  the 

functional  dependence  of  upon  n(r)  for  the  two  extreme 

limits,  it  is  probably  not  presumptuous  to  assume  that  for  a 

particular  type  of  aerosol  and  spectral  region  the  attenuation 

or  extinction  is  most  critically  dependent  upon  the  volume  of 

particulate  in  the  atmospheric  path  and  not  so  much  upon  the 

detailed  description  of  the  distribution  function.  For 

limited-visibility  atmospheric  conditions  we  have  generally 

found  a weaker  linear  correlation  between  extinction  at  or  near 

photopic  wavelengths  and  the  liquid  content.  As  expected,  there 

2/3 

is  a stronger  correlation  with  the  effective  area  density,  v 

For  example,  Mie  calculations  for  a set  of  measured  particle 

distributions  such  as  those  illustrated  in  Pig.  1 yield  a test- 

fit  relationship  8 ~ for  X = 0.55  unj. 

aer 

In  principle  and  perhaps  in  practice,  selection  of  a single 
moment  does  not  provide  a real  restriction  to  the  modeling  of 

y " " - ■ ■ 

The  r^  dependence  can  be  Justified  on  purely  physical  grounds 
by  noting  that  the  absorption  of  light  is  proportional  to  the 
cross  section  (area)  intercepted  times  the  path  length  tra- 
versed through  the  absorbing  mediiun. 
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infrared  propagation  in  aerosols.  For  example,  if  selection  of 
a single  moment  such  as  liquid  water  content  does  not  provide  an 
adequate  accounting  of  propagation  characteristics,  it  is  possi- 
ble to  carry  out  an  analysis  based  upon  two  or  more  moments. 

One  could  possibly  measure  or  monitor  the  moments  by  using  a 
multiple-fr^equency  light  detection  and  ranging  (LIDAR)  system. 

Using  Kle  calculations  for  58  different  particle  size 
distributions  representative  of  hazes,  fogs,  and  clouds  (Refs. 

3,  5,  b),  as  well  as  Independent  field  measurements  of  particle 
sizes  made  in  Grafenwohr,  FRG,  during  a 12-hour  period  (Ref.  2), 
we  have  establishsd  that  a strong  relationship  does  exist  be- 
tween the  total  volume  content  of  ^he  particulate  along  the 
transmission  path  and  the  aerosol  extinction  coefficient. 

Figure  3 illustrates  the  correlation  between  B at  10  pm  and 

aer 

liquid  content  and  the  subsequent  Independence  of  the  specifics 
of  n(r).  The  size  distributions  are  characterized  by  mode 
radii  ranging  over  approximately  four  orders  of  magnitude  up  to 
20  pm.  The  Grafenwohr  field  measurements  (•)  are  in  excellent 
agreement  with  the  theoretical  distributions  of  Refs.  3,  5 and 
6 (•). 

At  Grafenwohr  on  30-31  December  1975,  when  the  particle 
distributions  shown  in  Fig.  1 were  measured,  a simultaneous 
set  of  8-12  pm  extinction  measurements  was  carried  out  over  ar 
ll40-m  transmission  path  (Ref.  2).  The  straight  line  in  Fig, 

3 IS  the  regression  (p  = 0,88  for  12  data  points)  based  upon 
the  Grafenwohr  measured  size  distributions  (to  yield  the  mea- 
sured liquid  content)  and  measured  extinction.  The  relation- 
ship is  the  following: 


®8-12  «m  * ^39  » 


1.03 


where  v is  in  grams  per  cubic  meter.  This  is  in  very  good 
agreement  with  the  scaling  law  based  upon  Mie  calculations 
wltn  measured  or  assumed  distributions.  Since  our  Initial 
investigations  (Ref.  2)  suggesting  the  critical  role  of 


atmospheric  liquid  water  content  in  IR  propagation,  other 
workers  have  demonstrated  a similar  strong  correlation  betweo 
aerosol  extinction  and  liquid  water  content  for  maritime  er  '1 
ronments  as  well  (Ref.  7). 
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111.  CONCLUSIONS 

We  have  provided  a variety  of  experimental  and  theoretical 
data  to  demonstrate  that  to  a good  approximation  knowledge  of 
the  atmospheric  liquid  water  content  is  in  fact  sufficient  to 
determine  aerosol  extinction  of  the  far  Infrared  when  visibil- 
ity is  limited  by  water  hazes , fogs,  and  clouds.  Detailed 
knowledge  of  the  particle  size  distribution  is  unnecessary. 

Thus,  the  linear  dependence  of  upon  liquid  water  content 

suggests  that  a simple  measurement  of  a quantity  related  to  the 
volume  of  particulate  along  the  transmission  path  (l.e.,  grams 
per  cubic  meter)  could  provide  a direct  measure  of  the  IR  trans- 
mission characteristics. 
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